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Abstract 

We describe some new univariate spline quasi-interpolants on uniform parti- 
tions of bounded intervals. Then we give some applications to numerical analysis: 
integration, differentiation and approximation of zeros. 

AMS classification: 41A15, 65D07, 65D25, 65D32. 

1 Introduction 

Univariate spline quasi-interpolants (abbr. QIs) can be defined as operators of the form 

where {Bj,j G J} is the B-spline basis of some space of splines, say of degree d, on a 
bounded interval I = [a, b] endowed with some partition X n of / in n subintervals. We 
denote by 11^ the space of polynomials of total degree at most d. In general we impose 
that Q is exact on the space ILi, i.e. Qp = p for all p G ILi. Some authors impose 
further that Q is a projector on the space of splines itself (see e.g. [H],[I21)- As a 
consequence of this property, the approximation order is 0(h d+l ) on smooth functions, 
h being the maximum steplength of the partition X n . The coefficients fij are local linear 
functionals which are in general of one of the following types: 



(i) differential type : fij(f) is a linear combination of values of derivatives of /, of order 
at most d, at some point in supp(Bj) (see e.g. [TJ , [2] ) • The associated quasi-interpolant 
is called a differential quasi-interpolant (abbr. DQI). 

(ii) integral type : Hj(f) is a linear combination of weighted mean values of /, i.e. of 
quantities fwj where Wj can be, for example, a linear combination of B-splines (see 
e -g- [IB])- The associated quasi-interpolant is called an integral quasi-interpolant 
(abbr. iQI). 
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(iii) discrete type : f-ij(f) is a linear combination of discrete values of / at some points 
in the neighbourhood of supp(Bj) (see e.g. j2j, [U], [12], [El)- The associated quasi- 
interpolant is called a discrete quasi-interpolant (abbr. dQI). 

The main advantage of QIs is that they have a direct construction without solving any 
system of linear equations. Moreover, they are local, in the sense that the value of Qf(x) 
depends only on values of / in a neighbourhood of x. Finally, they have a rather small 
infinity norm, so they are nearly optimal approximants. 

In this paper, we only consider dQs, neither DQIs nor iQIs. We also restrict our study to 
splines defined on uniform partitions of I = [a, b}. Our aim is to give explicit formulas for 
dQIs of degrees 2 < d < 5 and some applications to three classical problems in numerical 
analysis: approximate integration and derivation and location of zeros of functions. The 
paper is organised as follows. In sections 2, we recall some facts about splines and 
quasi- interpolants. In section 3, we describe dQIs of degrees 2 < d < 5 and we give 
their infinity norms and their approximation orders. In sections 4 and 5, we give the 
associated quadrature formulas and some numerical examples. In section 6, we give 
the differentiation matrices for quadratic and cubic splines with numerical examples. 
Finally, in section 7, we show, on a simple polynomial example, how quadratic dQIs can 
be applied to the location of zeros of functions. 



2 Spline spaces on uniform partitions and dQIs 

For I = [a,b], we denote by Sd(I,X n ) the space of splines of degree d and class C d ~ l 
on the uniform partition X n = {xi = a + ih, < i < n} with meshlength h = 
A basis of this space is {Bj,j £ J}, with J = {1, 2, . . . , n + d}. With these notations, 
supp (Bj) = [xj-d-i, xj], and A/} = {xj-d, ■ ■ ■ , Xj-i} is the set of the d interior knots in 
the support of By As usual, we add multiple knots at the endpoints: a = x = x_i = 
. . . = and b = x n = x n+ i = . . . = x n+ d- We recall the represention of monomials in 
terms of symmetric functions of knots in J\fj |5].|19|. 

^6 { pBj{x), 6 { p = (^j symniriAfj), 0<r<d. 



In particular, the Greville points 



i=i 



are the coefficients of e\ = Ylj&j and the vertices of the control polygon of S = 
J2j£j c jBj are the control points V = {cj = (9j,Cj),j £ J}. The Schoenberg-Marsden 
operator is the simplest discrete quasi-interpolant which is exact on the space Iii: 



Sif = ^2f(e j )B j , s lP = P Vp£n x . 
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A discrete quasi-interpolant (abbr. dQI) of degree d is a spline operator of the form: 

whose coefficients are linear combinations of values of / on either the set 

T n ( for d even) or on the set X n (for d odd), where 

T n = {tj,j G J}, tj = + Xj-i), X n = {xj, 0<j< n}. 

Therefore, for d even, we set f(T n ) = {fj = f(tj), j G J}, and for d odd, we set 
f(X n ) = {fj = f{xj), < j < n}. Moreover, Q d is exact on ILj: 

QdP = P Vp G U d . 

For the construction of dQIs, i.e. for the determination of functional {/ij(f),j G J}, 
the exactness of Qd on 11^ amounts to solve a system of linear equations for interior 
B-splines and a finite number of specific linear systems for boundary B-splines. The 
determinants of these systems being Vandermonde determinants, there is existence and 
unicity of dQIs with the above assumptions (see also [3 j for more general cases). For the 
sake of completeness, we give below complete formulas for degrees 2 < d < 5. Moreover, 
we give exact values or upper bounds of ||Qd||oo and their approximation order on smooth 
functions. Actually, it is well known (see e.g. [H], chapter 5) that for any subinterval 
Ik — [xk-ii Xk], 1 < k < n, and for any function / 

11/ - Q d f\\oo,i k < (i + ||Q d ||oo)<W/, n d ) 

where the distance of / to polynomials is defined by 

^oo,/ fc (/,n d ) = inf{||/ -p\\oo,i k ,p e rii} 

Here, as usual, ||/ — p\\oo,i k = rnax x& i k \f(x) — p(x)\. Therefore, for / smooth enough, 
e.g. / G C d+1 (I), this implies that ||/ - QtfW^ = 0(h d+1 ). 

3 Discrete Quasi-Interpolants of degrees 2 < d < 5 
3.1 C 1 Quadratic dQI 

For the C 1 -quadratic dQI Q 2 f = S?=i A*i(/) Bji ^ ne coefficient functionals are easy to 
compute (details are given in [To].[T6]): 

»l(f) = h, &(f) = |(-2/r + 9/ 2 - / 3 ), /in+l(/) = |(-/n + 9/n +1 - 2/ n+2 ), 

A*n+a(/) = fn+i, and for 3 < j < n 

i*j(f) = k-fs-i+m-fj+i) 



3 



The exact value HQ2II00 — 1-4734 has been computed in [T6] . Therefore, for / e C 3 (I) 
for example, we have the following error estimates 

11/ - Q 2 /|U A < ^00,4 (/, n 2 ) for 1 < k < n =► || / - Q^IU = 0(h 3 ). 

3.2 C 2 Cubic dQI 

For the C 2 cubic dQI Q 3 f = Ysj=i /%(/) -^j> the coefficient functional are respectively: 
/i 2 (/) = i(7/ + 18/x - 9/ 2 + 2/3), /in +2 (/) = ^(2/„_ 3 - 9/„_ 2 + 18/ n _x + 7/ n ), 
= /o, Atn+3(/) = /n, and for 3 < j < n + 1 

A*i(/) = ^(-/i-s + 8/ i - 2 

As \n2\00 = \fi n +2\oo = 2 and \fMj\oo = | for 3 < j < n + 1, we obtain the upper 
bound || Q3 1| 00 < 2. It is possible to improve that result by writing the operator in the 
"quasi-Lagrange" form: 

71+3 

Qsf = E 

i=i 

where the fundamental functions are linear combinations of B-splines, e.g. for 4 < j < n, 
Lj = ~(— -Bj+3 + 8Bj +2 — Bj + %). It is well known that HQ3II00 is equal to the Chebyshev 
norm of the Lebesgue function: 

n+3 

A3 = £1^1- 

In each interval of the uniform partition, A3 is bounded above by the cubic polynomial 
whose Bernstein-Bezier (abbr. BB-) coefficients are sums of absolute values of BB- 
coefficients of fundamental functions. This allows to see that the maximum of A 3 is 
attained in the interval [a^a^] and we obtain: 

IIQ3II00 = IIAslloo « 1-631. 

From that we deduce for / G C 4 (I) for example, we have the following error estimates 

11/ - Qsf\\oo,l k < jj<W/, n 3 ) for 1 < k < n ||/ - Q3/IU = 0(h 4 ). 

3.3 Quartic dQI 

For the C 3 quartic dQI Q^f = Y^=x /■*?'(/) Bji the coefficient functionals are respectively: 

t*i(f) = fi, Vn+i(f) = fn+2, 

^a(/) = 105/1 + 32/2 - ge/s + 1B0/4 - 22l/5> 

HzU) = -45/1 + 28s/2 + 28s/3 _ 480^ + 288^, 
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^(/) - jkfi - + is/a - + 1^3/5, 

< f \ — 47 f 77 f I 251 f 97 J I 47 f 

rn+HJ j ~~ 315-'' 1 + 2 144^™+1 "r 144 •/« 240-> n ~ 1 1008-' n_2 ' 

/; r f _i£ / _L377f 1 J>lf _J>9_f I 7 j 

ri+2U J — ^45 Jn+2 T^gg./ri+l ^288-'" 2l 480'' n ~ 1 288^ n ^ 2 ' 

[J>n+3\f) = 195/1+2 + 32/n+l ~~ g^/™ 16o/ n_1 ~~ 224"/ ra_2 ' 

and for 5 < j < n 

Vjif) = YY^(/i-4 + fj) - 28g(/j- 3 + /i" 1 ) + 192/j- 2 

Let us give some details on the computation of functionals = 2, 3, 4 . As Hk( e r) = 
for < r < 4, we determine the five coefficients of the discrete functional 

fJ>k(f) = afe/i + Pkh + Ikh + S k U + Ckh 

as solutions of the three corresponding linear systems (2 < k < 4) of 5 x 5 linear 
equations 

t[a k + f 2 [3 k + tl lk + t\5 h + f 5 Ck = e k r) , < r < 4 

They have the same Vandermonde determinant V^ii, t 2 , ^3,^4,^5) 7^ since the t-s are 
distinct. Therefore they have unique solutions. The same technique is applied to the 
computation of other coefficient functionals. 

As I/X2I00 = l/Ai+sloo ~ 1.77, I/X3I00 — 1^+2(00 ~ 2.09, |/x 4 |oo = |AWi|oo ~ 2.88, \^j\oo ~ 
2.49 for 1 < j < 5, we can conclude that ||Q 4 ||oo < 2.88 and that for / e C 5 (7) for 

example, we have the following error estimates 

11/ - Q4/|U,4 < 4rfoo, 4 (/, n 4 ) for 1 < k < n 11/ - Q 4 /||oc = 0(h 5 ). 

3.4 Quintic dQI 

For the C 4 quintic dQI Q 5 / = Y^=i ^j(f) Bji the coefficient functionals are respectively: 
A*i — /o) AWs — /n, 

1^2 = |p/o + /l — J2 + 3/3 _ 7 l/ 4 + 2|/ 5 ^ 

A*3 = 200/0 + ^/i - I/2 + 10/3 - ^n/4 + ^k/5 

.. _ 41 f 1 43 f , 103 x 7 f , 13 f 13 r 
^ 4 ~ ~400/0 + go h + 120^ 2 ~~ 10^3 + 43/4 ~ 309/5 

( symmetric formulas for n + 2<j<n + 4), and for 5 < j < n: 

As l^loo = 1^+4(00 — 3.5, I/X3I00 = IAW3I00 ~ 3.92, I/X4I00 = IAW2I0C ~ 2.69, and 
\pLj\oo ~ 2.87 for 5 < j < n, we deduce that HQslloo < 3.92. Using a similar technique as 
for cubics, we find that HQ5II00 ~ 3.106. Therefore, for / G C 6 (I) for example, we have 
the following error estimates 

11/ - Q 5 f\\oo,i k < 4.5doo A (/, n 5 ) for 1 < k < n =}► || / - Q5/IU = 0(/i 6 ). 
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4 Application to numerical integration 

Newton-Cotes formulas are obtained by integrating interpolation polynomials (see e.g. 

[I! 0) In the same way, integrating spline quasi- interpolants give interesting quadra- 
ture formulas (abbr. QF) which are easily deduced from the above computations. We 
use the notations 

Af)= [ f, Uf)= f Qdf = Y,H(f) f Bj, E d (f)=i(f)-i d (f). 

J a J a j£J a 

As j^Bj = \{xj — Xj-d-i) and fij(f) are known explicitly, we can compute the follow- 
ing quadrature formulas. Moreover, as QIs give the best approximation order, we can 
conclude that E d (f) = 0(h d+1 ) for / G C d+1 (I), where h is the meshlength. Moreover, 
as for Newton-Cotes formulas, we get a higher approximation order for even degrees. 

(i) QF for quadratics 

MS) = JaQlf = h Si + h + M + l(h + /n+l) + I(/s + /„)] 

Error : for / G C\l), E 2 (f) = 1(f) -I 2 (f) = h*D*f(c) - ^h 5 D*f(c) 
This result is proved in [T6j . 

Error for Simpson: E*{f) = 1(f) - l*(f) = h A D A f(c) 

By comparing the two above errors, we see that the linear combination (extrapolation): 

X 2 (/) = l(32X 2 (/) + 23X*(/)) 
is such that 1(f) -l 2 (f) = 0(h 5 ). 

(ii) QF for cubics 

MS) = fiQ*f = hY£tfj+h [g(/0 + /») + l(fl + /n-l) + i(/2 + fn-2) + i(/ 3 + fn- 

Error: E 3 (f) = 1(f) — X 3 (/) = 0(h 4 ) for / G C A (I). Numerical experiments show that 
this formula is not as good as the preceding one. 

(iii) QF for quartics 

f b a QJ = MS) = h Eg Si + h [HOi + f n+2 ) + m ( / 3 + fn+l) + 1||( /3 + /b) ] 

+ h [fll^ 4 + + iH(/5 + /n-2)] 

Error: X(/) — X 4 (/) = 0(h 6 ) for / G C 6 (J). This is a remarkable formula, which can be 
compared to the Newton-Cotes formula of the same order. Numerical experiments show 
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that the error for the former QF has also the opposite sign of the error for the latter, as 
in the quadratic case. The proof will be given elsewhere. 

(iv) QF for quintics 

un = faQsf = iir; Id + * [K/o +/„) + ka + a-o + k/ 2 + /„_ 2 )] + 

Error: E 5 (f) = 1(f) — X 5 (/) = 0(h 6 ) for / G C 6 (J). Numerical experiments show that 
this formula is not as good as the preceding one. 



5 Numerical examples 

We compare numerical results on QF applied to the computation of 

1{f\) = J 1 - \q x 2 dx and ^(A) = J e~ x sin(5-Kx)dx. 
(i) QF/dQI degrees 2 and 3 

E 2 (f) = 1(f) - l 2 (f) = 0(h 4 ), E 3 (f) = 1(f) - l,(f) = 0(h*) for / G C\I) 
Simpson's QF E*(f) = 1(f) - l*(f) = 0(h 4 ) for / G C\I) 
Example 1: l(fi) 



n 


E* 




E2 




E3 




128 


0.73(- 


-9) 


-0.55(- 


-9) 


-0.44( 


-8) 


256 


0;45(- 


-10) 


-0.33(- 


-10) 


-0.26( 


-9) 


512 


0.28(- 


11) 


-0.21(- 


-11) 


-0.15(- 


-10) 


1024 


0.18(- 


12) 


-0.13(- 


-12) 


-0.95(- 


-12) 



Example 2: l(f 2 ) 



n 


E* 




E 2 




£■3 




128 


0.14(- 


-6) 


-0.11(- 


-6) 


-0.92(- 


-6) 


256 


0.90(- 


-8) 


-0.67(- 


-8) 


-0.52(- 


-7) 


512 


0.56(- 


-9) 


-0.41(- 


-9) 


-0.31(- 


-8) 


1024 


0.73(- 


-9) 


-0.52(- 


-9) 


-0.37(- 


-8) 



(ii) QF/dQI degree 4 
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E 4 (f)=l(f) -l 4 (f) = 0(h G ) for / G C\I). 

Newton-Cotes QF of degree 4: E* 4 (f) = 1(f) -11(f) = 0(h 6 ) for / G C 6 (J). 
Example 1: l(fi) 



Example 2: l(f 2 ) 



n 


E 4 


El 




128 


-0.83(-12) 


1.10(- 


12) 


256 


-0.12(-13) 


0.24(- 


13) 


512 


-0.18(-15) 


0.37(- 


15) 


1024 


-0.29(-17) 


0.59(- 


17) 


n 


E 4 


E\ 




128 


0.23(-7) 


-0.68(- 


-7) 


256 


0.44(-9) 


-1.04(- 


-9) 


512 


0.73(— 11) - 


-1.62(- 


11) 


1024 


0.12(-12) - 


-0.25(- 


12) 



(iii) QF/dQI degree 4: E 4 (f) = 1(f ) - l 4 (f) = 0(h 6 ) for / G C 6 (J). 
QF/dQI degree 5: E 5 (f) = 1(f) - l 5 (f) = 0(h e ) for / G C 6 (J), 
Example 1: l(fi) 



Example 2: l(f 2 ) 



n 


E 4 


E 5 




128 


-0.83(-12) 


0.95(- 


11) 


256 


-0.12(-13) 


0.14(- 


■12) 


512 


-0.18(-15) 


0.21(- 


■14) 


1024 


-0.29(-17) 


0.32(- 


■16) 


n 


E 4 


E 5 




128 


0.23(-7) 


-0.27(- 


-6) 


256 


0.44(-9) 


-0.50(- 


-8) 


512 


0.73(-ll) - 


-0.83(- 


■10) 


1024 


0.12(-12) • 


-0.13(- 


11) 



6 Application to numerical differentiation 

Differentiating interpolation polynomials leads to classical finite differences for the ap- 
proximate computation of derivatives. Therefore, it seems natural to approximate 
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derivatives of / by derivatives of Qdf as long as it is possible, i.e. up to the order 
d — 1. The general theory will be developed elsewhere. Here we only give results for 
the first derivative and d — 2, 3. We evaluate (Qdf)' = A t j(/)^ a ^ points T n for d 
even and at points X n for d odd. 

(i) Differentiation matrix for quadratics 

The derivation matrix T> 2 G R( n+2 ) x ( n + 2 ) is defined as follows: setting y G lR n+2 for the 
vector with components yj = f(tj),j G J and y' G IR n+2 for the vector with components 
Vj = (Q2f)'(tj),j G J, we simply write: 

y' = V 2 y 



v 2 = 



/- 


-8/3 


3 


-1/3 
















\ 




-7/6 


11/16 


13/24 


-1/16 


















1/6 


-3/4 


1/48 


5/8 


-1/16 


















1/16 


-5/8 





5/8 


-1/16 




















1/16 


-5/8 





5/8 


-1/16 


















1/16 


-5/8 


-1/48 


3/4 


-1/6 


















1/16 


-13/24 


-11/16 


7/6 


{ 



















1/3 


-3 


8/3 / 



(ii) Differentiation formula for cubics 

The derivation matrix G R( n+1 ) x(n+1 ) is defined as follows: setting y G IR n+1 for 
the vector with components yj = f(xj),0 < j < n and y' G IR n+1 for the vector with 
components y'j = (Q 3 f) f (xj), < j < n, we obtain: 









y' 


= V 3 y 










-11/6 


3 


-3/2 


1/3 













\ 


-1/3 


-1/2 


1 


-1/6 
















1/12 


-2/3 





2/3 


-1/12 
















1/12 


-2/3 





2/3 


-1/12 


















1/12 


-2/3 





2/3 


-1/12 
















1/12 


-2/3 





2/3 


-1/12 
















1/6 


-1 


1/2 


1/3 
















-1/3 


3/2 


-3 


11/6 J 



(iii) Some numerical results 

Again we use the two functions f\(x) = and f 2 {x) = e~ x sin(hnx) on the in- 

1 + IQx 2 

terval / = [—1, 1]. For p = 1, 2, we set e p = raax v£V \ f' p (v) — (Qdf P )'(v)\ where V — T n 
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(resp. V = X n ) for d — 2 (resp. d — 3) and £* = max„ e y |/p(i>) — 5/ p (i;)| where Sf p (v) 
is the classical centered approximation of f' p {v) of order 2 (with standard modifications 
at the endpoints). 

For quadratics, we obtain the following results 



n 


£\ 


F* 


£2 


F* 
£ 2 


64 


0.014009 


0.047853 


0.016143 


0.046317 


128 


0.003138 


0.012079 


0.003674 


0.011606 


256 


0.000767 


0.003036 


0.000872 


0.002902 


512 


0.000190 


0.000759 


0.000212 


0.00725 


1024 


0.0000475 


0.0001899 


0.000052 


0.000181 



We see that the orders are all 0(h 2 ). However, the errors for the derivatives of the 
quadratic QI (e\ and £3) are between 3 and 4 times less than the errors for the centered 
finite differences (e* and £3). 

For cuhics, we obtain the following results 



n 


£1 




F* 
fc l 




£2 




F* 
L 2 




64 


3.0(- 


-3) 


4.7(- 


-2) 


1.0(- 


-2) 


4.7(- 


-2) 


128 


2.0(- 


-4) 


1.2(- 


-2) 


1.4(- 


-3) 


1.2(- 


-2) 


256 


1.3(- 


-5) 


3.0(- 


-3) 


1.8(- 


-4) 


2.9(- 


-3) 


512 


8.0(- 


-7) 


7.6(- 


-4) 


2.4(- 


-5) 


7.2(- 


-4) 


1024 


5.0(- 


-8) 


1.9(- 


-4) 


3.0(- 


-6) 


1.8E( 


-4) 



Of course, e\ and £3 are both 0{h 2 ) and £1 and £3 are both at least 0(/i 3 ). However, 
for the function /1, a superconvergence phenomenon occurs because we have S\ = 0{h A ) 
instead of 0(h 3 ). We shall study this kind of results in a further paper. 



7 Approximating zeros of a function by those of a 
quadratic dQI 

Let / be a continuous function defined on I — [a, b]. In order to locate the zeros of / in 
this interval, we approximate / by its C 1 quadratic dQI g = Q 2 f and we compute the 
exact zeros of g: this is quite possible because g is piecewise quadratic. The complete 
study will be done elsewhere. Here we take a simple example: we want to approximate 
the zeros of the Legendre polynomial P%{x) = in the interval I — [—1,1]. The five zeros 
of P$ are respectively {±xi, ±#2, ^ x 3j ±^4}; with 

x l = .1834346425, x 2 = .5255324099, x 3 = .7966664774, x A = .9602898565. 
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The following array gives the errors Ek = Xk — Xk, 1 < k < A where Xk is the zero of g 
nearest to Xk- 

n e 1 e 2 e 3 e A 

16 .000543 .003784 .013753 -.007841 
32 -.000043 .000210 .000556 -.001017 
64 -.000013 -.000012 .000043 .000026 

Acknowledgements: the author thanks Professor Catterina Dagnino and the Depart- 
ment of Mathematics of the University of Turin for their kind invitation to deliver this 
seminar during his stay from January 12 to 20, 2005. 
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